The Angular Power Spectrum of the First- Year WMAP Data 

Reanalysed 



Pablo Fosalba 1 , Istvan Szapudi 



ABSTRACT 

We measure the angular power spectrum of the WMAP first-year temperature 
anisotropy maps. We use SpICE (Spatially Inhomogeneous Correlation Estima- 
tor) to estimate C/s for multipoles 1 = 2 — 900 from all possible cross- correlation 
channels. Except for the map-making stage, our measurements provide an in- 
dependent analysis of that by Hinshaw et al. (2003a). Despite the different 
methods used, there is virtually no difference between the two measurements for 
i 700 ; the highest £'s are still compatible within 1 — a errors. We use a novel 
intra-bin variance method to constrain Cg errors in a model independent way. 
Simulations show that our implementation of the technique is unbiased within 
1% for I ^ 100. When applied to WMAP data, the intra-bin variance estimator 
yields diagonal errors ~ 10% larger than those reported by the WMAP team for 
100 < i < 450. This translates into a 2.4 a detection of systematics since no dif- 
ference is expected between the SpICE and the WMAP team estimator window 
functions in this multipole range. With our measurement of the Cgs and errors, 
we get x 2 /d.o.f. = 1.042 for a best-fit ACDM model, which has a 14 % probabil- 
ity, whereas the WMAP team (Spergel et al. 2003) obtained x 2 /d.o.f . = 1.066, 
which has a 5 % probability. We assess the impact of our results on cosmological 
parameters using Markov Chain Monte Carlo simulations. From WMAP data 
alone, assuming spatially flat power law ACDM models, we obtain the reioniza- 
tion optical depth r = 0.145 ± 0.067, spectral index n s = 0.99 ± 0.04, Hubble 
constant h = 0.67 ± 0.05, baryon density Vt h h 2 = 0.0218 ± 0.0014, cold dark 
matter density Q c d m h 2 = 0.122 ± 0.018, and as = 0.92 ± 0.12, consistent with a 
reionization redshift z re = 16 ± 5 (68 % CL). 

Subject headings: cosmic microwave background — cosmology: theory — meth- 
ods: statistical 
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1. Introduction 

The Wilkinson Microwave Anisotropy Probe satellite (WMAP) has provided the clear- 
est view of the primordial universe to date. Its unprecedented high sensitivity and spatial 
resolution resulted in a unique set of cosmic microwave background (CMB) radiation maps 
with close to full sky coverage and uniformly high quality. As a result, fundamental cosmo- 
logical parameters can be constrained to the highest precision ever. Thorough analysis of this 
dataset (Bennett et al. 2003) yielded a cosmic variance limited measurement of the angular 
power spectrum, C/s, of the CMB temperature anisotropy for multipoles £ ^ 350 (Hinshaw 
et al. (2003a); hereafter H03). This confirmed and improved measurements from previous 
experiments (e.g., Miller et al. (1999); de Bernardis et al. (2000); Hanany et al. (2000); 
Halverson et al. (2002); Mason et al. (2003); Scott et al. (2003); Benoit et al. (2003)). The 
acoustic peak structure revealed by the WMAP temperature and polarization power spec- 
tra provided strong observational support to inflation and constrained viable cosmological 
scenarios to the domain of flat ACDM models and its close variants. 

Considering the importance of these results, our principal aim is to estimate the angular 
power spectrum in a completely independent way in the full range of multipoles probed by 
WMAP, 2 < £ < 900, and systematically compare results to H03. Our Ce estimation pipeline 
is based on SpICE 1 (Spatially Inhomogeneous Correlator Estimator; Szapudi et al. 2001a,b), 
a quadratic estimator based on correlation functions. SpICE performs edge corrections and 
heuristic minimum variance weighting in pixel 2 space to produce nearly optimal results. 
Our fast HEALPix 3 implementation of SpICE scales as (9(N 3//2 ) (N is number of pixels). 

2. Power Spectrum Estimation 

Our estimation methodology closely follows that of H03, but adapted to our technique: 

Step 1: We use the foreground cleaned intensity maps for the 3 highest frequency bands 
Q, V & W downloaded from the LAMBDA website 4 . Strong diffuse Galactic emission and 
resolved point sources are masked out using the KpO and Kp2 masks, that leaves 76.8% and 
85.0% of the sky useful for cosmological analyses, respectively. Monopole £ = and dipole 

1 htt p : / / www . ifa. hawaii . edu / cosmowave / 

2 The harmonic space alternative using pseudo C^s is MASTER (Hivon et al. 2002). 
3 http: / / www.eso.org/science/healpix/ 
4 http://lambda. gsfc.nasa.gov/ 
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£ = 1 terms are also removed from non-masked pixels. 

Step 2: Power-spectrum estimation is performed via SpICE: we compute the cross- 
correlations from 28 different pairs of channels constructed from the 8 "differencing as- 
semblies" (DAs) Ql through W4. Noise correlation among different channels is negligible, 
therefore our cross-power estimator is unbiased with respect to the noise (see e.g., H03). Like 
H03, we implement an heuristic ^-dependent pixel noise weighting scheme that minimizes 
errors: we use flat weights (mask weight only) for £ < 200, inverse pixel noise variance for 
£ > 450, and a transitional inverse rms noise weight in the intermediate range 200 < £ < 450. 

Step 3: A model for the power spectrum for unresolved extragalactic radio sources is 
subtracted from the cross-power spectrum of each channel. We implement the model given 
in §3.1 of H03. 

Step 4 : CV S from different channels are optimally combined using an inverse noise 
weighting, with DA sensitivities as described in the LAMBDA website. All channels are 
included, except for those in Q-band that are only used in the intermediate £-range. This 
helps minimizing galactic contamination at low £ and the window function cut-off at the 
highest multipoles. 

Step 5: Our quadratic estimator is defined in pixel space, where mask effects can be 
easily corrected for (cf. Szapudi et al. 2001a). The two point correlation function is then 
transformed into harmonic space via Gauss-Legendre quadrature to obtain the C/s de- 
convolved from the window function of the experiment. Symmetrized non-Gaussian beam 
transfer profiles (Page et al. 2003) and pixel window functions are corrected for in £-space. 

3. Principal Results 

Figure 1 shows the angular power spectrum of WMAP, AT 2 = £(£ + l)C^/27r, in /iK 2 
units, measured with SpICE. Upper panel shows the power spectrum for individual multi- 
poles, using Kp2 sky cut. Our measurement (red line) is in excellent agreement with H03 
(black line), multipole by multipole. In particular, for the quadrupole and octopole we find 
AT 2 ~ 135/iK 2 and AT 2 ~ 591/iK 2 , respectively (H03 get ~ 123/iK 2 and ~ 612/iK 2 ). For 
the highest ts we find slightly different amplitudes than H03, but consistent at the 1-a level. 

For the most part, we observe no systematic dependence of the measured CVs on the 
sky cut (see difference between red and blue lines in bottom panel of Figure 1). However, 
using KpO instead of Kp2 yields a 15% lower amplitude of the octopole £ = 3 and a 15 — 20% 
smaller amplitudes for the 3 highest band-powers centered at £ eS ~ 660, 750, 850. This effect 
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might be due to imperfect foreground removal and/or the intrinsic estimator variance due to 
finite volume and edge effects. We estimated the dispersion in a set of WMAP simulations 
with KpO & Kp2 sky cuts to be of the same order as the measured differences in the Cgs 
of the data. On the other hand, the cross-correlation amplitude between the clean WMAP 
maps and the best fit foreground templates is at the 5% and 10% level of the WMAP C/s 
for the lowest and highest ts, respectively. We thus conclude that sample variance due 
to sky coverage can account for most of the observed difference in the Q's, while residual 
foreground contamination is always subdominant. The low level of systematics in Kp2, and 
the increased statistical errors due to the decreased sky fraction left by KpO, motivate us to 
adopt Kp2 (as in H03) for the best estimate of the C/s. 



4. Error Estimation 

In order to estimate the covariance of our C/s, we generated MC simulations of the 
CMB sky and instrument noise for each of the 8 DA's (Ql through W4). We used the 
running index ACDM model that best fits a combination of WMAP, CBI & ACBAR data 
(denoted WMAPext in Spergel et al. (2003)). Maps were convolved with the symmetric 
(non- Gaussian) beam transfer function for each DA (Page et al. 2003). As for the noise 
simulations, we downloaded 100 sky maps per DA from the LAMBDA website. These 
simulate 1 full year of flight instrument noise and they include all known radiometric effects 
(Hinshaw et al. 2003b; Jarosik et al. 2003). Simulations were analyzed in exactly the same 
way as the data (see §2). All in all, we have constrained the errors from 1500 measurements in 
MC simulations (combining cross spectra using 100 MC's for each of the 6 highest frequency 
DA's VI through W4) for the multipole ranges £ < 200 and £ > 450, and 2800 measurements 
(combining cross-spectra using 100 MC's from each of the 8 DA's Ql through W4) for the 
intermediate Grange 200 < £ < 450. 

For multipoles £ < 350, errors in the WMAP power spectrum are dominated by cos- 
mic or sample variance (see H03) and the noise only contributes at the few percent level. 
Figure 2 displays the noise contribution to the relative errors at low multipoles, £ ^, 100. 
Correlated noise simulation results are displayed (smooth solid line) along with results from 
uncorrelated noise simulations (dashed line). The latter tends to underestimate errors by 
~ 1%. Alternatively the noise level can be estimated from the data rms dispersion among 
the WMAP channels used (oscillating solid line). These results are in excellent agreement 
with H03 (cf. lower panel in their Figure 4). 

At higher Ps pixel noise and systematic effects (e.g., beam and mode coupling, residual 
foregrounds) increasingly dominate the errors. MC methods assume detailed knowledge of 
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all such effects. To provide a model independent check of the errors, we introduce a novel 
technique that allows estimating errors directly from the data: the intra-bin variance (IBV) 
method. IBV estimates the variance of a given Cg from the rms dispersion in a bin B^ 
centered on I. The bin-width A£ is a matter of practical consideration, balancing variance 
and bias. More precisely, our estimator for a(Ci) reads 

a\C,) = ^— i Y,^Ce-^C i )f (1) 

where (ACg) = l/A£^, eB AQ/, ACg — Cg — Cf 1 , Cg is the mean of the measured Cgs 
over channels, and Cf 1 is our best guess for the data mean using a theoretical ACDM model. 
The latter is subtracted to decrease the bias due to the slope of the angular power spectrum. 
We used Cf 1 from the WMAP best-fit running index ACDM model (Spergel et al. 2003), 
although this is not critical: no baseline subtraction only biases at a few percent level. By 
construction, IBV should not be used to obtain errors with high resolution but to assess the 
overall level of errors in a range of ts, typically larger than Al. 

Figure 3 shows the ratio between the mean IBV rms dispersion to the usual MC rms 
dispersion, both estimated from ~ 3000 WMAP simulations of CMB signal and correlated 
noise. Narrow bins yield slightly biased (under-)estimates of the MC error at few percent 
level, possibly due to small mode-to-mode couplings. IBV method with Ai = 18 yields 
unbiased estimates of the error for WMAP simulations at the 1% level for £ £ 100. Doubling 
A£ introduces a slight high bias and significant edge effects for low ts that could be caused 
by the residual slope of the CgS. The unbiased bin-width, A£ = 18, with 15% variance is 
our choice for the WMAP error estimation. 

Figure 4 displays the WMAP data diagonal errorbars computed with the IBV method 
(spiky solid line) compared to the previously published diagonal errors (Verde et al. (2003), Hin- 
shaw et al. (2003b), Kogut et al. (2003)). The largest IBV errors appear to correlate well 
with the outliers of the data C/s with respect to the best-fit ACDM model (see Figure 3 in 
Spergel et al. (2003)), suggesting that our IBV estimator is closely related to a diagonal \ 2 
test. It is clear that the mean overall error is higher than originally estimated (otherwise the 
IBV curve would fluctuate around unity). The simplest and most conservative interpretation 
of our results yields a monotonic error increase with respect to the WMAP team diagonal 
errors of the form, a{Cg) IBV /cr{Cg) WM AP - 1-08 + 8.5 • 10~ 5 (£ - 100) for £ > 100 (straight 
solid line in Figure 4). This smooth prediction results from a least squares minimization 
to the IBV curve (large amplitude oscillating line in Figure 4). Note that for £ > 450, the 
error excess is consistent with the errors estimated from MC simulations with correlated 
instrument noise (see noisy line in Figure 4 growing from left to right). 

In the range 100 < £ < 450 the mean error level is incompatible with both MC simula- 
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tions that include correlated noise and the WMAP team published errors: given that there 
are approximately 16 independent A£ bins, with an intrinsic 15 % error each, and that the 
mean error excess is 9 % in this £-range, this amounts to a 2.4 a detection of the error excess. 
We have checked that using the C/s measured by Hinshaw et al. (2003a) yields comparable 
IBV errors in this multipole range. The excellent £-bj-£ agreement between the SpICE and 
WMAP team's measurement of the CVs (see Figure 1) indicates that both estimators win- 
dow functions are virtually identical in this regime, and thus the observed error excess points 
to systematics unaccounted for in the WMAP team analyses. For £ > 450 the interpretation 
is less clear as there are hints that both window functions might be slightly different (see 
lower panel in Figure 1). Such differences might arise in the practical implementation of the 
estimators. We also estimate a ~ 5% correlated noise contribution at £ < 100 (see Figure 2), 
that was neglected in previous likelihood analyses. A more robust assessment of errors is 
provided below using the full x 2 test, where off-diagonal terms are also taken into account 
following Eq.(15) in Verde et al. (2003). 



5. Discussion: Cosmological Parameters 

We investigate the implications of our measurements using a Bayesian analysis of cosmo- 
logical parameter estimation. We use CosmoMC 5 , a Markov Chain Monte Carlo (MCMC) 
implementation (Lewis and Bridle 2002) based on CAMB 6 (Lewis et al. (2000); see also CMB- 
FAST 7 , Seljak and Zaldarriaga (1996)). In order to allow direct comparison with Spergel 
et al. (2003), we focus on the simplest 6-parameter cosmological model consistent with the 
WMAP temperature and cross-polarization data. Following Verde et al. (2003), we assume 
a set of flat ACDM models with radiation, baryons, cold dark matter and cosmological con- 
stant. Primordial fluctuations are taken to be adiabatic and Gaussian with a power-law 
power spectrum. We use the physical dark matter Q cc i m h 2 and baryon VL^h 2 densities, the 
reionization optical depth r, the scalar spectral index n s , the normalized Hubble constant h, 
and the dark matter power spectrum normalization ag (Kosowsky et al. 2002). We estimate 
paramaters by combining 4 independent chains with 30000 accepted points each, and use the 
6 paramater covariance matrix as proposal density from precomputed runs. This yields an 
excellent convergence-mixing Gelman & Rubin statistic R — 1 ^ 0.02 for all cases studied. 

Table 1 summarizes our results. Imposing the prior r < 0.3, we find best fit values 



5 htt p : / / cosmologist . info / cosmomc 
6 http: / / camb.info 
7 http: / / cmbfast.org 
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Table 1. Best Fit Parameters for Power Law A CDM' 





SpICE C/'s 


WMAP CVs 




IBV Errors 6 


Standard Errors c 


T 


0.145 ±0.067 


0.151 ±0.069 


n s 


0.99 ±0.04 


0.99 ±0.04 


h 


0.67 ±0.05 


0.70 ±0.05 


tt b h 2 


0.0218 ±0.0014 


0.0234 ±0.0013 


h 2 

1 L cdm' L 


0.122 ±0.018 


0.123 ±0.017 




0.92 ±0.12 


0.92 ±0.11 


xtff/dof 


1398.8/1342 


1428.7/1342 



a WMAP Data Only. We impose a prior r < 
0.3. Table values are mean expectation values 
for the marginalized distribution and errors are 
the 68% (symmetrized) confidence intervals. 

b Parameters estimated with our MCMC's us- 
ing Ces measured with SpICE and IBV errors 
(see §4). 

c Same as b, but using C/s from H03 and errors 
from Verde et al. (2003); Hinshaw et al. (2003b); 
Kogut et al. (2003). 
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matching those of Spergel et al. (2003). In particular we obtain a x 2 /d.o.f. = 1.042 (i.e. 
it has a 14 % probability) for the best-fit model (see first column in Table 1), which is a 
slightly better fit to the data than that of Spergel et al. (2003), x 2 /d.o.f '. = 1.066 (i.e. 5 % 
probability). Our h and r are slightly lower but still consistent at the 1 — a level. This is 
more significant for our estimates of the C/s and errors (see first column in Table 1). In 
particular, our measurement flbh 2 = 0.0218 ± 0.0014 agrees with that from the latest BBN 
results tt b h 2 = 0.022 ± 0.002 (Cyburt et al. 2003; Vangioni-Flam et al. 2003; Cuoco et al. 
2003). We have checked that relaxing the r prior yields larger values of r = 0.19 ± 0.12 
(cf. Tegmark et al. 2004). Our main results (see first column in Table 1) are in excellent 
agreement with the best-fit values from WMAP+SDSS (Tegmark et al. 2004), and suggest a 
redshift of (abrupt) reionization z re = 16 ± 5 (68% CL). Data products and additional plots 
from this work can be found at http://www.ifa.hawaii.edu/cosmowave/wmap.html 
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sions, Olivier Dore, Hans K. Eriksen, Eiichiro Komatsu for useful comments, and Antony 
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the NASA Office of Space Science. Some of the results in this paper have been derived 
using HEALPix (Gorski et al. 1998). This research was supported by NASA through ATP 
NASA NAG5-12101 and AISR NAG5-11996, as well as by NSF grants AST02-06243 and 
ITR 1120201-128440. 
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Multipole 

Fig. 1. — WMAP angular power spectrum of the lst-year temperature anisotropy data. Up- 
per panel: Measurement for individual multipoles (red line) using the Kp2 sky cut. Results 
from by H03 are also shown (black line, nearly indistinguishable from red). The agreement is 
excellent for most £'s. Lower panel: Binned power spectrum for two different sky cuts, Kp2 
(red) & KpO (blue). We find excellent agreement with H03 (black line) for all multipoles 
I 700 and a slightly different amplitude for the highest £'s. 
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Fig. 2. — Errorbars estimated from uncorrelated (dash line) and correlated noise (solid 
smooth line) are shown. The spiky solid line shows the variance of the C/s among the data 
channels, what should depend on noise only. The Monte Carlo simulations with non-uniform 
noise appear to describe well the mean rms errors. 
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Fig. 3. — Solid lines show the ratio of intra-bin variance (IBV) to Monte Carlo errors 
measured in WMAP simulations. At i ~ 300 lines correspond to A£ = 6, 9, 18, 36 (bottom 
up). The third line (A£ = 18) is unbiased for i > 100. Upper dotted lines (in reverse order 
to solid lines) correspond to the relative rms error on the estimated IBV errors (e.g., in the 
convention used, 1.15 means 15% error). 
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Fig. 4. — Ratio of measured errors vs. published WMAP errors (Verde et al. 2003; Hinshaw 
et al. 2003b; Kogut et al. 2003). The noisy blue line and smooth red line through the latter 
show the usual MC dispersion and the average IBV error respectively, for the same set of 
simulations. The large amplitude oscillating black line shows the IBV estimator applied to 
the WMAP data. The IBV method detects 8-15% larger errors than the published ones. 
A conservative IBV error estimation is depicted by the solid straight line obtained from a 
linear least square fit. 



